Non-Crystalline Structures of Ultra-Thin Unsupported Nanowires 
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Computer simulations suggest that ultrathin metal wires should develop exotic, non-crystalline 
stable atomic structures, once their diameter decreases below a critical size of the order of a few 
atomic spacings. The new structures, whose details depend upon the material and the wire thickness, 
may be dominated by icosahedral packings. Helical, spiral-structured wires with multi-atom pitches 
are also predicted. The phenomenon, analogous to the appearance of icosahedral and other non- 
crystalline shapes in small clusters, can be rationalized in terms of surface energy anisotropy and 
optimal packing. 

PACS numbers: 61.46. +w, 61.43. Bn, 81.05.Ys 
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The regular crystalline packing of atoms, which usu- 
ally minimizes energy in a bulk solid, is known not to last 
indefinitely as physical size is decreased. For very small 
elemental clusters, for example, one expects a structural 
switch from bulk-like to some new arrangement at some 
critical radius [|IJ. Such is the case, for example, of small 
(diameter ~ 20-40 A) clusters of Ag Q and other metals 
||, and of rare gas clusters J|j5], which were observed 
to abandon single crystal structures in favor of an icosa- 
hedral shape below a critical size. Cluster calculations 
based on the Lennard- Jones (LJ) potential confirm 
the stability of the icosahedral shape, implying a radical 
deviation from the bulk fee morphology, for N < 4000. 
The driving force causing this morphological change at 
small size can generally be traced back to two very dis- 
tinct sources: 1) electronic magic sizes, stabilized by 
filling of shells, as in atoms and nuclei; 2) competition 
between optimal internal packing and minimal surface 
energy, the latter dominating for sufficiently small size. 
Here we shall address exclusively phenomena caused by 
the second mechanism, which is important in a large class 
of cases. For example, very small icosahedral clusters, 
rather common in nature, are favored by their ability 
to expose an exceptional 100% of (111) faces, lowest in 
energy. 

We address the same structural question for ultra-thin 
nanowires. A priori, one could anticipate a weaker ten- 
dency to depart from the bulk structure than in clus- 
ters, because of the smaller relative surface/bulk ratio 
in the wires (2/R against 3/R, where R is the radius). 
However, we present here detailed optimization results, 
suggesting that new non-crystalline nanowire structures 
(which we propose, for convenience, to name weird wires) 
should in fact be rather readily realized in nature. As it 



turns out, the tendency for the wire structure to switch 
from crystalline to weird appears to be both general, and 
stronger than anticipated. Moreover, there seems to be 
a variety of new structures, changing continuously and 
unpredictably from size to size, and from a substance to 
another. For our initial analysis we have chosen to study 
elemental (unsupported) metal wires, whose fabrication 
should be within reach, thanks to the metal's combined 
ductility and strength. 

We considered so far Al and Pb as test cases. Both 
were assumed to be described by classical, empirical 
many-body interatomic potentials, of the "glue" type [fr|. 
The glue parameters were previously optimized to model 
the properties of metallic Al m and Pb |)| respectively. 
To be sure, such potentials are not meant to provide a 
quantitative description for metal wires of nearly atomic 
thickness, whose electronic structure may and will dif- 
fer profoundly from that of the bulk metal. The struc- 
tures thus generated should therefore be considered at 
this stage of mostly qualitative significance, while future 
ab initio studies will be called for to achieve quantitative 
accuracy, and to include electronic shell effects, expected 
to be very important at small radii. Nonetheless, the ex- 
istence of stable weird wire structures, suggested by our 
energy optimizations with classical potentials, represents 
in our view a true, general, and parameter-independent 
phenomenon. For instance, it is found that even a hy- 
pothetical argon wire (as described by a Lennard- Jones 
potential) turns weird when sufficiently thin JToj] . We 
also note that the simplicity of the empirical potentials 
for metals has been at this initial stage quite instrumen- 
tal, as it allowed a broad and thorough search for optimal 
structures, which would be prohibitive, had we adopted 
from the outset the more accurate ab initio approach. 



1 



For structural optimization, we adopted a molecular- 
dynamics based "simulated annealing" methodology, us- 
ing techniques previously employed to investigate pre- 
melting of thin wires We generally started with 

fee (llO)-oriented TV-atom wires [in a few cases we also 
used (100) wires, with similar outcomes], with periodic 
boundary conditions (PBC) along the wire axis direction 
z. We first relaxed the atomic positions, to optimize the 
fee wire structure. Starting from this structure, we simu- 
lated thermal annealing cycles, with T rising in steps up 
to a maximum (typically about 0.75 T m , where T m is the 
bulk melting point) and then decreasing back to T = 0. 
The wire length was clamped by the PBCs, and its value 
kept constant during annealing, thus preventing surface- 
free-energy-driven contraction into a drop, which thermal 
diffusion would otherwise inevitably drive. Once back at 
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T = 0, however, the length was finally allowed to adjust, 
and the wire structure further relaxed, before inspection 
of its energy and structure. The resulting annealed wire 
energy was almost invariably better than that of the fully 
relaxed original fee structure. Its shape generally exhib- 
ited large irregularities and a ^-dependent thickness, oc- 
casionally with sections showing a new regular structure, 
which initially triggered our curiosity. Simple annealing 
was generally unable to generate a fully regular wire. In a 
generality of cases where the spontaneous regular sections 
could be identified by visual inspection, new full-length 
entirely regular wires were produced with just that struc- 
ture. In the new artificial wire, the total length and atom 
number were chosen such as to guarantee a proper match- 
ing through PBCs. After further relaxation we obtained 
the final wires, collectively depicted in fig. |l|. 
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FIG. 1. Total energy per atom E vs. inverse wire radius 1/R for the relaxed structures obtained by optimization for Al 
(a) and Pb (b) wires. A selection of morphologies is shown, fee wires are represented by open squares, "weird" wires by full 
circles. Very thin wires which do not belong to either class have been marked with open triangles. Solid lines represent a fit 
to fee wires using eq. (jjj), and dashed lines a fit to weird wires using (Q). Weird structures become favored for R < R c . Inset: 
structure of the helical Pb wire B7: complete wire (n = 1, 2, 3), with outer shell removed (n = 1, 2), and inner strand (n = 1). 
Note how the outer shell exhibits a nearly-square atomic structure, while that of the second shell is nearly-triangular. Also 
note the different helical pitches of these two shells. 



The wire radius R is defined by S = irR 2 , where S is 
the area projected on the xy plane by all the atoms (con- 
sidered as spheres of diameter d, where d is the nearest- 
neighbor distance in the bulk crystal) included in a por- 
tion of the wire of length 2<i-\/2/3, corresponding to two 
(111) interplanar spacings. 

Crystalline wires clearly prevail for R larger than a 



critical value R c of the order of 3d, depending on the 
system. For R < R c , weird structures make their ap- 
pearance. Their energies fall systematically below the ex- 
trapolated fee line, consistent with simulations, showing 
an irreversible and exothermic spontaneous restructur- 
ing of the initially fee wire. The structures found for Al 
and Pb are found to bear certain resemblances, but are 
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certainly not identical, although we have not attempted 
a complete search. A (simplified) visual illustration of 
some stable weird wire morphologies is offered in Fig. [|. 
These structures, though noncrystalline, are very regu- 
lar. The thicker ones, about four atomic radii in diame- 
ter, are three-shell wires (A15, B6, B7), where a central 
string is surrounded by two successive coaxial cylindrical 
shells, mutually related by some kind of epitaxy. Several 
structures, including three-shell (A15, B6) and two-shell 
(A6, B4) cases and several composites (A8-A14), display 
pentagonal motives. 

For better characterization, we have calculated angu- 
lar correlation functions for all our structures. These are 
defined by considering all the angles formed by bonds ij 
and ik, where i runs on all the atoms closest to the wire 
axis, and j and k on their neighbors. The cutoff distance 
used to define the neighbors is 3.3 A for Al and 3.5 A for 
Pb. Results for a selection of cases are shown in fig. [|. 
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FIG. 2. Angular correlation functions, restricted to corner 
atoms located at or near the wire axis. Al wires are on the left 
panel, and Pb wires on the right panel. Dashed lines corre- 
spond to perfect icosahedral angles (63.4°, 116.6° and 180°). 
A16 and B8 represent a nearly perfect fee packing, and A6 is 
a nice example of icosahedral order. 

Wires A16 and B8 exhibit peaks at 60, 90, 120 and 180 
degrees typical of regular fee structure. In contrast, the 
pattern of the pentagonal wire A6 has only three peaks 
near 63, 117 and 180 degrees, the signature of icosahe- 
dral packing. Strong icosahedral features can be also 
observed on All and A15. A local icosahedral geome- 
try has recently been found in a Na neck studied by first 
principles molecular dynamics [l3| . Occasionally we find 
structures like A10 exhibiting a broad distribution, indi- 
cating a marginal case. In the thinnest wires (probably 
less meaningful than the others), we also observe trian- 
gular and tetrahedral packings. Pb wires also exhibit 
pentagonal motifs (B4, B6), but angles show no icosahe- 
dral packing: this different behavior is discussed below. 



Some wires, such as Al, A5, A9, A12, B7, are truly 
weird, possessing chiral, or helical, arrangements, with 
long periods. In these cases we also checked for stabil- 
ity against a discrete change of pitch, although we could 
not generally rule out the possibility of a continuously 
changing (incommensurate) pitch. 

The detailed internal structure of a beautiful weird 
wire, the three-shell helical B7 of Pb, is illustrated in 
Fig. |l|b (inset). The center (n — 1) consists of a single, 
nearly straight atomic chain, or strand. The surrounding 
shells, n = 2, 3, can be seen as being made up of m n iden- 
tical strands, each helically wound with a linear pitch A„ 
in a cylinder of radius r n . Choosing conventionally the 
strand direction closest to the wire axis z, we follow each 
strand and determine the number M n of periodically re- 
peated cells and the corresponding number of turns N n 
required for that strand to connect back to itself. Hence 
A„ = M n c/N n , where c = 60.83 A is our cell length. For 
wire B7 we find m 2 = 7, r 2 = 3.3 A, M 2 = 7, N 2 = 2, 
A 2 = (7/2)c, and m 3 = 11, r 3 = 6.0 A, M 3 = 11, N 3 = 7, 
A3 = (ll/7)c. Strikingly, the two consecutive shells are 
not simply commensurate, which implies in this case an 
overall coincidence length of (7 x ll)c, or 4684 A! 

How should we understand these new morphologies? 
The basic explanation is clearly related to surface energy. 
In the weird wires, contrary to fee wires, the "bulk" pack- 
ing is not so good, but the surface packing is excellent. 
Unlike crystalline wires (which must by necessity possess 
different crystallographic facets), the weird wires in fact 
display a single, optimal surface structure. This feature 
is also a characteristic of, e.g., icosahedral clusters, and 
represents the winning ingredient at very small radii. 

Quantitatively, we can write the crystalline wire energy 
per atom as 



£ c (i?) = -£ coh + ^7- 
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where -E C oh is the bulk cohesive energy, il c = nR 2 L/N 
(L is the wire length) is the average atomic volume, 7 
is the surface energy (averaged on the different facets), 
and ft is the average edge energy per unit length (there 
are 8 edges in our geometry). 7 can be written as 
7 = 7111 [am + ai o/ioo + aiio/uo] where 7m is the 
close-packed surface energy, a%jh is the fraction of ex- 
posed area relative to (ijk) facets, and fijk — 7yfc/7m- 
Moreover, 8ft = 4/xi + 4^2 where fix and fi 2 are the edge 
energies relative to the (lll)-(lOO) and (111)— (110) edge 
respectively. The weird wire energy, on the other hand, 
is much simpler: 
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reflecting a) a single smooth closed-packed surface and b) 
no edges. This surface gain is countered by bad packing, 
causing a bulk energy increase A and an average atomic 
volume fl w (defined as for the crystalline wires) slightly 
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larger than £l c . The crystalline- weird transition will take 
place at 



1 

R c 
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where R c0 — ^Sflcft/irA, and S — 7 — (Sl 11J /ri c )'7 111 
is an effective surface energy difference taking into ac- 
count the different atomic volumes, that can be either 
positive or negative, depending on the material. From 
separate calculations, we know a priori all parameters, 
except for A and ft. The geometries of the sections of 
fee wires are adequately represented by am — 0.55, 
aioo = 0.25 and a no = 0.20. For Al (Pb), fi c = 16.1 
(27.9) A 3 , Q w /n c = 1.10 (1.00), 7111 = 54.3 (37.5) 
meV/A 2 , / 100 = 1.08 (1.01), / uo = 1.19 (1.11), 7 = 57.5 
(38.5) meV/A 2 . 

Figure [I] shows the resulting plots for eq. ([!]) , obtained 
with ft as a fit parameter, and of (§) with A as the sole 
parameter. The accord with energies obtained by simu- 
lations is quite good, yielding A = 37 (24) meV, ft — 35 
(25) meV/A, and R c = 5.3 (9.8) A. These values of A 
are in good agreement with the T = energy difference 
between the glassy metal (obtained by simulated quench- 
ing of the liquid) and the fee solid, namely 51 (29) meV. 
The weird wire interior is therefore close to a glass. The 
edge energies are also reasonable, since they should be 
close to half the step value, or about jd/A, which is 39 
(33) meV/A. We conclude that eq. (||) should have pre- 
dictive power for other materials, provided the glass, the 
surface and the edge energies can be estimated. 

This physics also explains the diversity of behaviour 
between the two metals Al and Pb, both fee in bulk. 
The surface anisotropy of Pb is much smaller than that 
of Al. All stable Al wires display an outer shell which is 
triangular — like a curved (111) surface — to a very good 
degree. In Pb, conversely, weird wires tend to possess a 
roughly square outer layer, which can be explained with 
the exceptionally low 7100 0]. For instance, the pen- 
tagonal (non-icosahedral) structure B4 with (100)-like 
lateral faces is preferred over A6 (icosahedral), similar 
but with (lll)-like lateral faces. In the helical wire B7, 
the presence of a relative tilt angle between third-shell 
and second-shell strands, 63 — 62 = tan _1 (27rr3/A3) — 
tan~ 1 (27rr2/A2) — 16.0°, appears to realize a good ap- 
proximate coincidence of the third (outer) shell atoms 
with the hollow sites of the second shell, which in turn 
has a roughly triangular packing (inset). This provides 
an amusing case of what might be considered "curved 
surface epitaxy" between two otherwise incompatible 2D 
lattices. 

Finally, we compare wires with clusters. The non- 
fec wire structures we just found for Pb are plentiful. 
Conversely, the same potential is known to stabilize only 
fee crystalline clusters and to destabilize icosahedra, in 



agreement with experimental indications |TJ], down to 
the lowest size (13 atoms) ||. Hence, the tendency of 
wires to abandon the bulk fee structure is much stronger 
than could have been expected, in particular stronger 
than in clusters. The reason is most likely related to 
the well-known fact that icosahedral shapes require (111) 
faces to be stretched in order to form a space filling 
structure ||. Owing to the high tensile surface stress 
of the metal, this circumstance can disfavor icosahedral 
clusters against fee structures, when the surface energy 
anisotropy is small, as in Pb. On the other hand, this 
negative factor is absent in wires, which are open, and 
thus can fully adjust along z. A second difference is that 
icosahedral clusters have edges, while most weird wires 
have none. In both respects, therefore, wires are not nec- 
essarily just a two-dimensional version of clusters, and 
new phenomena can arise. 

Could such unsupported metal nanowires be fabri- 
cated, and their possibly weird structure eventually be 
detected? Wires of over a thousand Angstroms length 
can be pulled by an STM tip |l6) . They can be expected 
to possess, at least in some sections, the ultra-small radii 
which we have addressed here. Weird structures could 
be sought in field ionization or transmission electron mi- 
croscopy (TEM) images. Also, the electronic structure 
of a wire is in principle a strong function of its atomic 
structure. Future study might reveal a measurable im- 
print of its weird shape, if present. More work in these 
directions is currently planned. 
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